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Abstract 

We propose a family of boundary terms for the action of a causal set with a 
spacelike boundary. We show that in the continuum limit one recovers the 
Gibbons-Hawking-York boundary term in the mean. We also calculate the 
continuum limit of the mean causal set action for an Alexandrov interval in 
flat spacetime. We find that it is equal to the volume of the codimension-2 
intersection of the two light-cone boundaries of the interval. 
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1 Introduction 

One approach to constructing a quantum dynamics for the causal set approach to 
quantum gravity [1] is to discover a discrete counterpart of the gravitational action, 5[C] 
that can furnish the weight, e lS ^ c \ of each causal set, C , in the gravitational sum over 
histories. A start in this direction has been made with a proposal for scalar curvature 
estimators for causal sets of dimension d [2-4], Summing such a scalar curvature 
estimator over all elements of a causal set (causet for short) gives a natural proposal for 
a causet analogue of the Einstein-Hilbert action, a proposal that remains to be studied 
in depth. The purpose of this paper is to initiate an investigation of gravitational 
boundary terms for the action of causets. This is likely to be important as, in the 
continuum, it is well known that the Einstein-Hilbert action, Seh, is not the full story 
in the presence of spacetime boundaries. Indeed, the gravitational action must include 
a boundary term Sghy, the Gibbons-Hawking-York (GHY) boundary term, in order 
to yield a well-defined variational principle when the metric is fixed on the boundary 
of spacetime [5, 6]. If the classical limit of quantum gravity is to arise from the path 
integral in the expected way, such a term in the action will be essential when boundaries 
are present. Whilst we do not yet know how to fix boundary conditions for causets 
in general, it is likely to be useful to have an analogue of the GHY boundary term 
for any causal set which is well-approximated by a manifold with a boundary. In this 
paper we initiate the study of boundary terms for causal sets by proposing an analogue 
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in the case of spacelike boundaries and investigating the above mentioned causal set 
action for a causal interval with null boundaries. First we consider causal sets which are 
well approximated by (M,g), a d-dimensional, causal, Lorentzian spacetime with finite 
volume which admits a closed, compact spacelike submanifold, E, such that the causal 
past and future sets, M ± := J^E), satisfy M + D M~ = E. Then E is a component of 
the future (past) spacelike boundary for M~ (M + ) and the GHY term for E, considered 
as a boundary of M + or M - , is given by 

Sghy^M*] =Tj^ Jj^xVhK , ( 1 ) 

where K is the trace of the extrinsic curvature K p n a of E, l p = ( 87 rG)^ 
is the rationalised Planck length and we are working in units where H — 1. Here we 
take n p to be the future-pointing timelike unit covector normal to E. We work with a 
mostly plus convention for the metric so is past-pointing. 

We recall that the integral in (1) is equal to the normal derivative of the volume of 
E along the unit normal vector field, nG 

f d d ~ 1 xVhK=-^~ f d d ^ 1 xVh, (2) 

Jr, on 

where this is the rate of change of the volume backwards in time, as is past-pointing. 
This observation suggests a natural candidate for an analogue of the GHY boundary 
term for E for a causet that can be faithfully embedded in M. Spacetime volume 
corresponds to cardinality in a causet. Hence the spatial volume gradient corresponds 
intuitively to the difference between the number of causet elements that are future 
nearest neighbours of E and the number of past nearest neighbours. This intuition 
turns out to be a good guide and we will identify a family of causal set boundary terms 
based on it. The family of causet boundary terms we find corresponds to the different 
ways to define a discrete derivative that tend to the same limit in the continuum. We 
will define the causet functions and state the claims about their continuum limits in 
section 2. The calculations justifying the claims appear in section 3. In section 4 we 
investigate the proposed causal set action for causal intervals in flat spacetime and show 
that its mean takes the form, in the continuum limit, of a boundary contribution from 
the codimension 2 “joint” of the interval’s boundary. 
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2 The Boundary Terms 


Given a finite causet ( T >, ^), we define an u F k element” to be an element of T> with 
exactly k elements - not equal to itself - to its future in the order 1 and define F k [D ] 
to be the number of F k elements in V. Similarly a u P k element” is an element with 
exactly k elements to its past in the order and P k [D ] is the number of P k elements in 
V. For example, Fq[D] (P 0 [^]) is the number of maximal (minimal) elements of V. 

Given a finite causet, (C, z<), with two subcausets, C + and C~ we introduce the 
following family of causal set “boundary terms” (CBT): 


S^ BT {C,C-x*-,p,q] := (i/y 


d—2 






PmPm 


[c-] + E^.n 


where the constant ad is given by 


d(d + 1) / s d . 2 y 
(d + 2 ) Vd(d-i); 


( 3 ) 


( 4 ) 


S d = (d + 1 ) 7 T~ 2 ~/r (4id + l) is the volume of the unit d-sphere, l is a length and l p 
the Planck length defined previously, p and q denote finite strings of real numbers 
(po, • • • ,Pmi ■ ■ ■) an d (qo ,..., q n ,.. .) respectively. We will prove that the strings must 
satisfy the following conditions: 


r(4 + m ) r (4 +n) 

Y^Pm U ^ ^ = 0 , 

ml 


\ ^ r ( 2 d + m) \ - 

yPm -j-V Qn 

ml 


n\ 

r(! + n) 


= 1. 


n\ 


( 5 ) 

( 6 ) 


We call (3) a boundary term but in general, when C + and C~ are arbitrary subcausets 
of C, it will have no physical significance. 

Let (M, g) be a d-dimensional spacetime with finite volume and spacelike, closed, 
compact hypersurface S as described in the introduction. Given such a spacetime (M, g ) 
and sets M ± := J ± (E), S^ BT defines a family of random variables in the following way. 
The Poisson process of sprinkling points into M with density p = F d generates a random 
causet (C, together with subcausets C ± which consist of those elements sprinkled into 
M ± . The functions P k and F k acting on the random causets C + and C~ respectively 

1 Other ways to express this are, “has exactly k descendants” and “precedes exactly k elements 

excluding itself”. 
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are random variables P*. and F k . These random variables can be substituted into (3) 
to give the family of random variables S^ T : 


s cbt l M ,z,p;p,q] : = ( l / l p) d 2 a d 




(7) 


We claim that- in the limit of infinite density the expectation value, in the sprinkling 
process, of S^ BT tends to the continuum GHY boundary term of the surface £: 


lim 

i-s-0 



1 


Id- 2 
L P 


d 


d-1 


X 


Vh K 


Sghy 


[e, m~] , 


( 8 ) 


where (•) denotes the mean over sprinklings. We will prove this in the next section. 

One can see that at least two non-zero entries in p and q together are necessary in 
order to satisfy (5) and (6) and if exactly two entries are non-zero they will be uniquely 
fixed, but if more than two entries are non-zero this uniqueness is lost. This accords 
with the continuum boundary term being a first derivative. The freedom of choice in p 
and q is the freedom to discretise a derivative in many ways but the difference of two 
nearby values is sufficient. 

We introduce special notation for the simplest member of the family: 

sf [C,£T,C + ] := (F„ [C-] - P 0 [C+]) . (9) 

This is proportional to the difference in the numbers of minimal elements of C + and 
maximal elements of C~, and thus corresponds to the intuitive idea described in the 
introduction. An illustrative sketch of the idea is shown in Figure 1. This case is the 
easiest to investigate computationally, and we shall use its random variable counterpart, 

[M, E,p], later when we study the fluctuations of the discrete boundary terms 
numerically. 

There are two special subfamilies of boundary terms, one defined by p = 0 and the 
other by q = 0. In the former (latter) case, this corresponds to defining a boundary 
term for the past (future) boundary of C + ( C~ ) using only data from C + (C~) itself. 
The simplest cases of these boundary terms are 

:= m) d - 2 Ah (Po [C + }-dP, [C+]) , (10) 

. ( 11 ) 

These give rise to random variables [M, Y>, p\ and [M, E, p\ via sprinkling at 

density p = l~ d as before. 
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Figure 1: An illustration of a sprinkling into a spacetime partitioned by a spacelike 
hypersurface. Black points correspond to causal set elements and links (irreducible 
causal relations) between elements are shown as thin black lines. The maximal (F 0 ) 
elements in C~ and the minimal (P 0 ) elements in C + have been highlighted with white 
hlling. The shaded areas illustrate the regions whose volumes V k and Vj are needed in 
the proof of Section 3. In this sketch Pq = 11 and Pq = 5 (with time flowing upwards). 


2.1 The Surface Volume Family 


We also propose a family of causet functions that will give the volume of a spacclikc 
hypersurface in the appropriate context: 

A^[C,C~,C + -,p,q\ ■= {l/l p ) d ~ 1 b d [Y J Pm F m [C-]+J2 ( ^ P n [ C+ ] ) > ( 12 ) 


where 


and p and q now satisfy 


bd = d 


S r , 


d—2 


d(d — 1) 


r (h + m) r (1 + n) 

X>m i +Eg" ■ = 1 • 


n\ 


(13) 


(14) 


We see that only one non-zero entry is necessary to give an expression for the discrete 
surface volume. Once again, for (M, g), E and p = l ~ d , we can define a family of random 
variables, 


A^ [M, E, p:p, q] := (l/l p ) d ~ 2 b d f ^p m F m + ^ q n P„ j 

\ m n J 


(15) 


We claim that, in the limit of infinite density, the expectation value of A (</) in the 
sprinkling process tends to the spatial volume of the surface E: 

1 


&o ( A(d) [ M,S,P; ^^l) = pW / d<1 lx ^h- 


(16) 
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In proving (8) we will establish the results necessary to prove this claim as well. 

One can define functions for the volumes of future and past boundaries respectively 
as the two simplest members of the family: 


Af[C-\ ■= Wl P ) d -' P/y F„[C-] , 

(17) 

A^[C + ] := (W J_I jttL a[C + ] . 

(18) 


3 Calculations 

3.1 (F m ) and (P m ) for Poisson Sprinklings 

In order to establish (8) and (16) we will find the behaviour of the mean values of F m 
and P m as p tends to infinity (/ —> 0). The result will follow almost immediately. 

For any given realisation of the sprinkling, the probability that a sprinkled point 
p e M~ below the surface E is an F k element of M~ is given by the probability that k 
points of the sprinkling he in the region J + (p) D J _ (E). The Poisson process assigns a 
probability, 

P (k points in J + (p ) fl J~(E)) = e~ pVy( ' p ' > , (19) 

to this event, where V T (p) := vol (J + (p) D J~( E)) is the spacetime volume of the region 
J + (p ) D J“(E). The probability of sprinkling an element into an infinitesimal d -volume 
dV p at p is pdVp where p = l~ d is the sprinkling density, and so the total expected 
number of F k elements below E is 

(Ffc) = p f dV p ^ML e ~ P v,{p) ( 20 ) 

Similarly the expected number of P k elements above E is 

(P fc ) = P f dV p (PV !,[ P)) e, (21) 

7j+(s) fc! 

where V k (p) := vol(J + (E) fl J~(p)). 

Let = (t,x) be “synchronous” or Gaussian Normal Coordinates (GNCs) adapted 
to E such that in a neighbourhood C/s of E the line element is 

ds 2 = —dt 2 + hij(t,x.)dx l dx j . (22) 

In these coordinates the surface E corresponds to t — 0. Each point p in the neighbour¬ 
hood lies on a unique timelike geodesic whose tangent vector at E is equal to — id L on E, 
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and the t coordinate of p is equal to the proper time from E to p along that geodesic. 
The spacetime (M, g) restricted to this neighbourhood of E is globally hyperbolic and 
E is a Cauchy surface within it. 

The integrals (20) and (21) seem intractable as they stand, since the integration is 
over the entire causal past/future of the surface. However, since E is closed and compact 
and M + and M~ are of finite volume, we can always find a subneighbourhood of Us such 
that the contribution to the integrals from the complement of that subneighbourhood 
tends to zero exponentially quickly as p —>• oo. Let e > 0 be small enough such that for 
all p G E and \t\ < e, (t, x(p)) are the GNCs of a point in LG. Define Us(e) := {q G 
Us : |f(g)| < e} and consider the integral in (21) restricted to W := J + (E) \ Us{s)'- 


dV n 


<J>VM) k „- P v> 


(p) 


’w 


k\ 


< 


-pVl 


dVr, 


I w 


{pVMt 

k\ 


(23) 


where ||e -py *|| is the uniform norm over the integration region W. Since V A (p) increases 
with t along the geodesics from E, and { x p E Us ■ t = e} is homeomorphic to E and so 
is closed and compact, V A (p) achieves its minimum value V min > 0 in W at some point 
with t — e. Then ||e _pVi || = e~ pVmin and so the integral (23) falls off exponentially fast 
as p —> oo. Similarly for (20). 

Thus, so long as p is large enough, we make only an exponentially small error by 
cutting off the integration ranges in (20) and (21) at t — ±£ with £ as small as we need 
in order to be able to expand in powers of t. Expanding the determinant of the metric 
around t — 0, the integrals we want to evaluate are 


(F k ) = P I 

f d d ~ x x I 

f dt hh ( 

J 

s J 

V 

< p *> = P ) 

f d d ~ x x \ 

f dth} ( 

J 

s J 

0 V 

where both h 

= det (hij 

) and h 


ft3 (l + + 0(f 2 )^| + 


2 h 


2 h 


k\ 


k\ 


(24) 


dh 

dt 


that vanish exponentially fast in the limit p —> oo”. 


3.2 Volumes 

Let us first calculate V A (p ), the volume of region X p := J~{p) flM + , for p in Us{s) flM + . 
Let po be the point on E where the unique timelike geodesic through p whose tangent 
is normal to E intersects E. Let the values of p’s GNCs be x p = (£ p , x p ), then 
po has GNCs £q := .x' /J '(po) = (0, x p ). We choose £ small enough that, for every 
p G Us(s) flM + , there exists a Riemann normal neighbourhood centred on p 0 containing 








the region J~(p ) D M + . We choose Riemann Normal Coordinates (RNCs) centered at 
Po, y^ = (y°,y) = (t, y), such that the GNC time coordinate of p equals the RNC time 
coordinate of p: t p = t p =: T. 

The relationship between RNCs y M and GNCs x v is, to second order, 

t = AK*" - <) + - x^x" - 4) + 0(0c - l0 ) 3 ) . (25) 

The constant matrix A ^ obeys 

A^A v v ri w = g pu (po) , (26) 

and the metric and Christoffel symbols in RNCs are flat at po: 

9w(Po) = Pw > 
r^(Po) = o. 

The inverse coordinate transformation is 

x^^ + //+0(|/ 3 ) , (28) 

where A C is the inverse matrix of A^, i.e. A^A% = 5^ and A^A^ = 5^. There is no 
0(y 2 ) term in (28) due to the fact that TCRpo) = 0. The components of A ^ satisfy 
A° 0 = 1, A°i = 0, A^A^Sy = hij{po ) and = A\A^h i:j (p 0 ). 

As we saw above, the limit we are considering, p —y oo, allows us to take £ to be 
arbitrarily small so the limit can be thought of as driving each relevant p —y p 0 on E 
along the geodesic normal to E. That makes the region X P) whose volume we need, 
tend to a truncated solid, nearly flat cone with apex p and a base on E defined by a 
quadratic form in the three spatial RNCs around Pq. The leading contribution to the 
volume, 

V k {p)= [ d d y \/—g{y) , (29) 

Jx P 

is therefore the volume, T d , of the flat cone of height T with a flat base on surface 
t — 0. Corrections to this are higher order in T and come from three sources: (i) 

V~9(y) + !, (ii) the null geodesics down from p to E forming the top boundary, 

T p := dJ~(p) fl M + , of X p are not straight, and (in) the base, B p := ED J~(p) is not a 
flat disc. The first two corrections are due to the curvature of M and the third comes 
from the extrinsic curvature of E. 

The correction from (Hi) is found by taking spacetime to be flat, so that RNCs 
are the usual Cartesian coordinates centred at po and T p is the top boundary of the 
flat cone, satisfying (v 1 ) 2 = (T — t) 2 and t £ [0,T], The base B p in GNCs lies in 
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surface t — 0, so we can use (25) to find the equation for the surface in RNCs. This 


gives 

i = ^r° J .(p 0 )(^ - xi)(x j - x 3 0 ) + 0{{x - x 0 f) . (30) 

The linear part on the right hand side of (25) vanishes, since A 0 (x 11 — Xq) = t (which 
follows from A° t = 0 and A° 0 = 1) and t = 0 on the bottom surface. Using the inverse 
RNC relation (28), the equation for B p in RNCs is 

t = )r»,(p„)V i V # V + 0(y 3 ) . (31) 

Let us rewrite this equation in spherical polar coordinates, i.e. define r := a/ 
and the usual angular coordinates 0i, (ftd -2 hi terms of the spatial coordinates y 1 = 
rcos(0i),... ,?/ d_1 = rsin(0i) ■ • • sin(0 d _ 3 ) sin(0 d _ 2 ). Then 

* = l (r“,(Po)V,V,^) r 2 + 0(y 3 ) = i/(x r , <p)r 2 + 0(y 3 ) , (32) 


where (p stands collectively for all the angular coordinates </>i, pd- 2 - The function 

/(x p , <p) depends on x p since T° • and A\ depend on p 0 . 

With the boundaries of the integration region in hand, we can now write down the 
integral explicitly in spherical coordinates: 

P r f'f'max (</>) r-r+T 

/ d d y = dn d _ 2 r d ~ 2 dr / dt + 0(T d+2 ), (33) 

JX p in flat space JS d ~ 2 J 0 J ^ f (x p ,(/>)r 2 

where r max (<p) is the value of the radial coordinate for which B p intersects T p at an 
angle 0, as shown in Figure 2. Equating the time coordinates of T p and B p gives 

^/(x p , 4>)r max 2 (4>) = -r max (0) + T . (34) 

We solve this for r max (<fi) and take the positive solution. The solution can be expanded 
in T and is simply r max = T + 0(T 2 ), with angular dependent terms contributing at 
0(T 2 ). The 0(T 2 ) term will contribute at 0(T d+2 ) in the volume integral. Substituting 
fmax = T into (33) allows us to evaluate the integral (33), which equals 


Sd-2 j,d 

d(d — 1) 


(l - 2pyy r «(j>o)V.' 4 y' ,T ) + 0 ( T “ +2 ) ■ 


(35) 


where the 5 T] comes from the fact that cross terms (i ^ j) vanish under the angular 
integration. The defining relations for A\ can be rearranged to give A l t A 3 -S^ = h l i(p 0 ), 
and in GNCs the extrinsic curvature on the surface is given by 


K = 


1 h 

2 h ' 


(36) 
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Figure 2: A 3-dimensional representation of the region X p in RNCs. The top, T p , of 
dX v can be approximated as a fiat cone, and the base, B p , intersects T p at a radial 
coordinate, r max which will in general be a function of the angles (p (in 3 dimensions 
there is one angle 0). This function is found by projecting down from the intersection 
to the t plane. 

Substituting this into (35) we obtain 

W^Ty Td 0 + 2PTT) A '<°- ^ T ) + 0(Ti+2 > • < 37 > 

and we see the first contribution is the volume of the flat cone with flat base as expected, 
and the first correction is of order T d+1 . 

The corrections (i) and (ii) come from the non-flatness of the metric. The deter¬ 
minant y/—g can be expanded in RNCs and the deviation of T p from straight lines 
considered. The curvature contribution to the volume of a small, approximately flat 
causal interval - or Alexandrov neighbourhood - of these effects has been calculated 
[7-9] and the same arguments show that the corrections (i) and (n) in our case are of the 
same order, 0(T d+2 ), which means they are suppressed with respect to the correction 
derived above. This is to be expected on dimensional grounds as extrinsic curvature 
has dimensions of inverse length whereas Riemann curvature has dimensions of inverse 
length squared. We will see that ( T d+2 ) corrections do not contribute to the boundary 
term in the limit. 

So we have 

^ r ' x > = wk) Td 0 + 2PTT) A '<°' X > T ) + 0{T “ +2 > ■ < 38 > 

The same calculation for p e U-z{e) fl M~ gives 

= i^r- rf“ t 1 - WTTy R ' (0 ' X)T ) + 0(TM) ■ (39) 
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where T > 0. 


3.3 The Mean of S 


(d) 

CBT 


Focussing on (P*,), we use (38) in (24) to find 

p k+i A k r - re 


<p*> = 


pk+2 

k\ 


d^x h ~2 / dt [( t dk + (kB - K) t dk+l ) + O ( t dk+2 )] e~ pAt 

Jo 

[ d d -'xhi [ dt +0(t llM+2 )]e-' ,A ‘ d + ... , 


(40) 


where K := A'(0,x) and we have defined 

A : = 

B : = 


S d . 


d- 2 


d(d~ 1) ’ 
d 


(41) 


-K . 


2 (d + 1) 

To find out how this behaves in the limit of p —» oo, it suffices to find the behaviour of 
the following integral: 

ff [ dt t q e~ pMd , (42) 

Jo 

where p, q G M. We make the substitution z = pAt d to put (42) into the form of an 
incomplete gamma function. 

A-m 


d 




dpAe d 


dz z( 2 ^)“ 1 e^ 


As p —> oo we have 


rpAe d 


i ! (¥K- = rfiilW..., 


(43) 


(44) 


where, as before, + ... denotes terms that tend to zero exponentially fast. 
Then the limiting behaviour of (P*,) is 


<p*> =p 1 ~ i (fa) -1 r - 


cP-'x Vh 


P 1 - 1 * (ad )' 1 r ^ | k ^ f d^xVhK + O^p 1 - i) , 


/ F \ (h \-i r 6 + fc ) 

(hfc)-p d (M --- 


(45) 


+ P 1- ^ M' 1 [ d d ~ 1 x VhK + O , 

* j Yj 
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and we have given the behaviour of (F*.) as well for completeness. 

The results for the boundary terms and the surface volume follow almost immedi¬ 
ately: 




hm {S { cbt) = ^ a d ( ^ Pm (F m ) + ^ q n (P„) 

” m n 


= lim 

p—t OO 


P* a d I V- „ T (l + m ) . V" rr r (5 + U ) 

l^bA^” 1 ml 


d d ~ l x Vh 


n\ 


l P~ 2 V 
P \ rn 

u 


T r (i +m ) - y>Pi±d) [ d ‘-i x vKk + o(p-i) 

ml n\ / ly 

m n / J ~ 

d d ~ l x VhK , 


(46) 


using the conditions (5) and (6) for p and q. Also 


lim ( A (,/) ) = lim 

p—yoo 


p 


4-1 


p ^00 ] 


-bd ( ^ Pm (F m) + ^ <?n (P n) 


= lim 

p —^00 


Id -1 




T 


(47) 


1 x A , 


using (14). 

The form of (45) suggests that (F*.) and (P*,) can be written as a Laurent series in 
the discreteness length, l, starting at l 1 ~ d . If the coefficients of the higher order terms 
in l are proportional to higher order normal derivatives of the surface volume then one 
might be able to find causal set analogues for second and higher normal derivatives of 
the surface volume. 

It is interesting to note that from the limiting behaviour of (P*.) or (F*.) we can 
find the dimension, d, by taking the ratio of either (P 0 ) and (Pi), or (F 0 ) and (Fi). 
The limiting behaviour of the latter ratio, expressed in terms of the discreteness length 
l, is found to be 


(Fo) = d _ farffi S^'xVKk 

(Fi> «,,I (4- 1 ) ps-'xVh 


(48) 


In the limit of l —4 0 one gets the dimension exactly. The fraction involving the two 
integrals is simply the average value of the extrinsic curvature across E. The case for 
the ratio of (P 0 ) and (Pi) is the same as (48) but with a positive sign after d. 
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3.4 Finite p and Fluctuations 

To decide under what circumstances the causal set boundary term, evaluated on a single 
causal set sprinkled into M, is close to the continuum GHY boundary term of E, it 
is necessary to know both the size of the fluctuations about the mean and when that 
mean is close to its limiting value. 

To take the second point first, the mean is close to its limiting value when the 
next order term in the expansions performed in the previous section can be ignored. 
Firstly, p must be large enough that an £ > 0 exists such that the expansions in GNCs 
are valid in a neighbourhood U^(e), and such that pV min 1 so e~ pVrnin <C 1, and 
the integral over the region outside U^(e) is negligible. V min ~ e d , and so e l. The 
expansions in equations (24), (38) and (39) are valid if at each point of E there exist 
RNC compatible with the GNC such that ICe <C 1 and Tie 2 <C He, where /C and 7 1 
stand for any component of the extrinsic curvature of E and spacetime curvature of 
M, respectively, evaluated on E. The resulting conditions are Til 2 <C HI <C 1 which is 
just what one would expect. If these conditions did not hold then a discrete manifold 
with discreteness on scale l could not be expected to encode the geometry of E and M 
around E. 

We now turn to the fluctuations or standard deviation, cr[Sp^ T ] = Var[S^g T ]^, of 
the causal set boundary term around the mean. A heuristic argument gives an estimate 
of the dependence of fluctuations on p = l~ d . In any spacetime region of fixed volume 
V the number of causal set elements in a sprinkling is a random variable, N, with 
mean (N) = pV and s.d. \J (N). Consider the simplest boundary term S^. The 
volume of a region corresponding to a thickening of the hypersurface E by one unit 
of the discreteness scale l (e.g. by Lie dragging the surface along its normal by an 
amount l ) is approximately vol(E)/ = vol(E)p _ a. Since Fo and Po are random variables 
that count nearest neighbours of E we may therefore expect their mean values to scale 

d —1 d —1 

like pvol(E)/ = vol(E)p“ oc (N) d , and indeed this agrees with the leading order 
behaviour of (45). This suggests that P 0 and F 0 will be subject to fluctuations of order 

d— 1 d— 1 

(N) 2d = [pV)^r in the limit of large p. Moreover F 0 and Po are independent and so 
cr[S^ T ] should behave like p^rp^r = phsr . Hence for d = 2 these fluctuations should 
grow like pi as p —> oo, for d = 3 they should be constant, and for d > 3 they should 
be damped. 

We tested this with simulations in the simplest case of flat spacetime and flat 
surface E. We took a sample of 100, density p = l~ d , sprinklings of a h-cube [0, \} d 
in d-dimensional Minkowski space with hypersurface E : t — 1/2, and evaluated the 
sample mean and (corrected) sample standard deviation of The expectation value 
of S{n is exactly zero due to the symmetry of the situation. 


14 



6 


m 

o 

CO 


CM 

O) 

o 



10 15 

l_og 2 N 


20 


Figure 3: A plot of the standard deviation in samples of 100 of for a flat (K = 0) 
surface bisecting a d = 2, 3 and 4-dimensional unit cube in Minkowski space for different 
values of N = p. Black dots, blue triangles and red squares correspond to the simulation 
results in d — 2,3 and 4 dimensions, respectively. The corresponding black, blue and 
red lines have gradients |, 0 and and best-fit intercepts of order 1. 


Figure 3 shows the results for d = 2,3,4 spacetime dimensions, with (N) = p 
ranging up to 2 20 . Each data point represents the sample standard deviation for a 
sample of 100. The solid lines have been obtained by fitting an arbitrary constant 

3-d 

multiplier in the scaling law predicted by the argument above, 7(d) x (N) 2d , to the 
data. The best fit values are all of order 1: 7(2) = 0.80, 7(3) = 0.97, and 7(4) = 1.07. 
The data are evidence for the scaling predicted by the heuristic argument. The sample 
means (not shown) for different p are consistent with zero within the standard error. 
Simulations for the boundary term (which is proportional to dF\ — Fo) show the 
same dimension dependent scaling behaviour for the standard deviation, though in this 
case the heuristic argument is complicated by the fact that the random variables of 
which the boundary term is a sum are not independent. 2 

2 While the heuristic argument predicts a scaling of the mean and standard deviations consistent 
with the data, a closer look at the samples we generated for and P*, for k = 0,1 suggests that their 
distributions deviate from a Poisson distribution: they are “underdispersed”, i.e. their s.d. grows like 
the square root of the mean but is related to it by a constant of proportionality less than 1. We have 
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4 The Causal Set Action for a Flat Alexandrov In¬ 
terval 


Now that we have a family of analogue GHY boundary terms for causal sets in hand we 
can consider if such terms need to be included in any putative action for causal sets. 
In particular we can ask whether boundary terms need to be added to the recently 
proposed Benincasa-Dowker-Glaser (BDG) causal set actions [2-4], Before that question 
can be answered, it is necessary to determine whether the BDG actions already contain 
any boundary contributions. 

The BDG action S^ DG [C] of a finite causal set C is 

\sf D a [C] = - ad (l/l r ) d - 2 U [C] + cfAqcd , (49) 

' i= 1 ' 


where N t [C] is the number of (i + l)-element inclusive order intervals in C, N[C] is the 
cardinality of the causal set, and l/l p is the ratio of a fundamental length to the Planck 
length 3 . The constants are 


and 


Old = 


r (i + i) 

rh+1)' 

d + 1 


: 2 / d d odd 


2/d 


2/d 


d even , 


_ 2+T (1 + |) d 

0d ■ ) r (j + 2)r(j + i) m 

r (i)r W J 


d odd 


d even , 



d odd 
d even , 


(50) 


(51) 


(52) 


where q = 2 1 iSd- 2 /{d(d — 1)) (recall that Sd is the volume of the unit h-sphere). The 
coefficients of the terms N t [C\ in the sum are 


- '1 £ (Ilf + 1 ) + i) 
hi U ^d + DGh + i) 

G, r(j(fc + i) + 2) 

hi U jr(f + 2)r(^ + i) 


d odd 


d even . 


(53) 


begun to investigate this further and hope to return to a more careful study of the distributions of 

these random variables in a future note. 

3 We reintroduce h in this section. 
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We note here that these coefficients can be expressed more compactly as generalised 
hypergeometric functions of type {q + 1 ,q}: 

C\ d) = q+ iF q ({ai,...,a q ,i-l},{bi,...,b q }\l) , (54) 

with q = a* = 4±A and b t — ^ for d odd, and q = |, a, = d+2 ^ +2 and h = ^ for d 
even. 

As in Section 2, given a causal Lorentzian spacetime (M,g), the sprinkling process 
at density p = l~ d turns this function of causal sets into a random variable SfocMp], 
the “random discrete action” of ( M,g ). A requirement for the causal set action to be 
physically interesting is that its mean should tend to the continuum action of (M, g) as 
p —» oo. The question at hand is whether in this limit it includes boundary contributions 
in addition to the Einstein-Hilbert term. 

We will explore this question by calculating the mean of the d- dimensional BDG 
action for causal sets sprinkled into causal intervals or “Alexandrov intervals” in d- 
dimensional flat spacetime. Since the Einstein-Hilbert contribution is expected to 
be zero, this will teach us something about what boundary contributions, if any, are 
included in the BDG action. The boundary of an Alexandrov interval consists of a 
past and a future null cone which intersect at a codimension-2 joint of topology S d ~ 2 
(see Figure 4). While the GHY term is defined in the continuum for both spacelike 
and timelike boundaries, and contributions from their codimension-2 intersections or 
joints have also been worked out [10], null boundaries and their intersections, on the 
other hand, are little discussed and there is no consensus on whether the variational 
principle for GR can be made well-defined on a region with null boundaries (see, however, 
[ 11 , 12 ]). 

As was shown in [13], when d = 2 the continuum limit of the expectation value 
of the discrete random action of an Alexandrov interval of arbitrary size is equal to 2. 
While this might suggest topological invariance, we will now show that it is a part of a 
more general result for d > 2 and has a geometrical origin. Namely, it corresponds to the 
volume of the joint of the Alexandrov interval, which in flat spacetime is a (d — 2)-sphere 
and independent of the interval size only in d = 2. 

Consider an Alexandrov interval, I(p, q), of proper height r between two points p 
and q in d- dimensional Minkowski spacetime. Its boundary consists of the two null cones 
from p and q which intersect at the joint, j( d ~ 2 l ■= dJ + (p) D dJ~(q), a codimension-2 
sphere of radius r/2. The joint has volume vol = (t/ 2) d ~ 2 Sd-2- The interval 

itself has volume vol (I(p, q)) = 2(Sd_ 2 /(<i(d — l)))(r/2) d . For the sprinkling process at 
density p = l~ d , the mean, N := (N), of the number of causal set elements sprinkled 
into I(p, q) is N = pvol(/(p, q)). In what follows we take the continuum limit p —> oo 
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Figure 4: The Alexandrov interval I(p,q). The boundary consists of the null sections 
B ± and the spatial sphere S d ~ 2 at their joint. 

while keeping r fixed. The mean of the random discrete action of this flat region should 
give, in the limit of large p, contributions from the boundary only. 

In [14] a closed form expression was obtained for the mean value of the number of 
(■i + l)-element inclusive intervals contained in an Alexandrov interval in d-dimensional 
flat spacetime: 

, ( «\ r(ci ) 2 ;v‘ +2 y. (—JV)*r (fc + j + i)r +1) r (gj±n +1) 

' * ' r W ho r (fc + i + 3) r (fc +1) r + d) r ( d(fc+ 2 m) + d) ’ 

where i > 1. Importantly, this power series can be expressed more compactly in terms 
of a generalised hypergeometric function of type {d, d} as shown in [14], and is therefore 
convergent for all N. All the power series in N that appear subsequently in this section 
are therefore also convergent. We now use this to evaluate (S^ DG ) in an Alexandrov 
interval in flat spacetimes of different dimensions. 

We begin with the simplest case of d — 2, where 

\ {s ( bdg) = 2 (n - 2 (n< 2) ) + 4 (n< 2) ) - 2 (n« 2 >)]) . (56) 
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Using (55) gives a power series expansion in N with coefficients — , i — 1,... oo, so 
that 

\ (Sbdg) = 2 (1 - e ~ N ) , (57) 

which agrees with the result in [13]. In anticipation of the results for higher d we note 
that the volume of the zero sphere at the joint, vol(,jU 0 ' ) ) = S 0 = 2, so that 

IhU ( S «7c) = ™l( J |0 >) . (58) 

This is in agreement with the result obtained for a 2-dimensional flat causal interval 
[13], 

Next, substituting (55) into the d — 3 averaged BDG action, 



gives a power series expansion in N with coefficients 

, (L\ H )" 1 8 

a \i P ) x a (3i + i)(3i-i) ’ 

where i — 1,... oo. Rearranging indices we find a closed form for the action: 


1 

h 


a(3) 

>bdg 



([-1 + 2^2 




(59) 


(60) 


(61) 


where 2 F 2 is a generalised hypergeometric function of type {2,2}. This can be re¬ 
expressed more simply as 



where y(s,x) = t s ~ 1 e~ t dt is a lower incomplete Gamma function. The large N 
behaviour is thus dominated by the last term in the above expression. Using 7 (s,x) = 
T(s) — T(s, x), where the upper incomplete Gamma function T(s,a;) ~ x s ~ 1 e~ x in 
the asymptotic limit, the dominant term in (62) simplifies to —Aa^lN 1 ^'^(2/3)/l p = 
vol (J^)/l p . Hence 

lirn \(S% G ) = 7 L vol(J (1) ) • ( 63 ) 

00 tl L p 

For d — 4 
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Excluding the first term, this is a power series in N with coefficients 


fiY v ( 3 !) 2 (-l) M (i-l)( 2 i- 3 )! 
‘\lj 3 i!( 2 i + l)! 


(65) 


where now i — 2,..., oo. Using this, Mathematica yields the closed form expression 


1 

h 


3 ( 4 ) 

’SDG 



— 3 (7 — e~ N + r( 0 , N) + ln(iV))^ , 


( 66 ) 


where 7 is the Euler-Mascheroni constant, and Erf is the error function. Since 
Erf(-s/]V) goes to f in the asymptotic limit, the dominant contribution to the 
above expression comes from the second term, —3aJ 2 V^ iV/Zp, which simplifies to 
2y/6n~Nl 2 /l 2 = vol (J 2 )/l 2 . Thus, again 

4A i ( S «U) = 5 vol( ' Jt2)) • (67) 

We now turn to the case of general d. We begin by writing the (averaged) sum in 
(49) as a power series in N: 

rid 00 

E cf ’ <N,> = E A f " i+1 ■ < 68 ) 

»=1 j =1 

After a rearrangement and redefinition of indices we find that 



2 (-it r(f(j -1) +1) r(fj +1) 
U +1)! r (5O' -1) + d) r (jjj + <i) 


D +2 


E(-d‘ 


1=1 


c. 


id) 


(69) 


where d = 2D for d even and d = 2D +1 for d odd. While (68) can be directly evaluated 
by Mathematica for small values of d — 2,..., 5, it is greatly assisted by the following 
simplifications to the for higher d. 

We begin by evaluating the sum in (69). We first use Mathematica to evaluate it 
for d — 2,..., 20 which then suggests the general form 


D +2 


J 


e(-du:;'^= 

i =1 V 


(-1)” ((2^ + 1 Yf - 1) (3 - (2D + l)j/2)p_ l 
4r(2 + D) 

(—l) D Dj(2 + 2D){1 — Pj) D -\ 

2Y{2 + D) 


d odd 

d even , 
(70) 
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where (a)& is the Pochhammer symbol. Inserting this into (69) we use Mathematica 
to evaluate it for d = 2,... ,20. After some manipulations this suggests the general 
expression 



r (d ) 2 (—1 )^ +1 

r (il) + 1 )) r (K 2 + j)) r ( 2 + t> 



(71) 


where 



>r(2 + dj) 

2l+dj p (d-k) 

r (1 + If) f (2 + if) 

r (i) 


d odd 

(72) 

d even . 


Taking our cue from the behaviour of (S^ DG ) for d 
we will consider the ratio 


2, 3,4 in the N —> oo limit, 


id —2 / a(d) \ 
L p \ a BDGJ 

h vo \(j( d ~ 2 ')) 


_ £_d _ 

d(d- l)r (1 + |) 


N+—J2 c i d) ( N *> 

ad *=i 


n-d 


(73) 


where £ ( i — 1 for d odd and 2 for d even. Finally inserting (71) into (73) Mathematica 
gives for d — 2,..., 16 

h ( S bdg) = jpVo\(j( d ~ 2) ) . (74) 

This is the main result of this section and can be interpreted as saying that, in the 
continuum limit, the mean of the random discrete action of a causal diamond is a 
pure boundary term coming only from the volume of the codimension-2 joint. One 
might speculate that the BDG action for the Alexandrov interval contains all the GHY 
contributions which implies that the null boundary GHY term vanishes identically. 
Interestingly, this concides with the claim in [11] that the GHY term for an Alexandrov 
interval is given only by the volume of the spatial joint. The result we have obtained 
is for flat spacetime and it would be interesting to see how the presence of curvature 
affects it by repeating this calculation in RNCs to the lowest order corrections. 

Finally, while efforts have been made to find a closed form expression of i(S ^ DG ) 
for arbitrary d this has proved difficult, even in the asymptotic limit. As we now show, 
the most obvious approach of using the asymptotic form of the (N-^)s is insufficient 
for this purpose. In the large N limit [14] 



r(l + »)r(d) 

(l-i) (I +1)« 


N [ 


2-3 


+ 0(N a f(N)) , 


(75) 
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where a — 1 for d — 3,4 and 2 — | for d > 4. and /(IV) = In IV for d = 4 and 1 otherwise. 
For d = 2 

^Nf } ^ = 7VlnIV + 0(IV) . (76) 

Since this dominates the leading order contribution of N^sr to \(S^ DG ) for all d, it is 
clear that this contribution must vanish. Inserting (75) and (76) into the BDG action 
confirms that this is indeed the case for d — 2,... 16. In fact, the next to leading order 
terms in (75) and (76) also do not have the requisite N~ dependence, and are dominant 
in comparison. Hence their contribution too should vanish, but we do not have an 
explicit expression for their coefficients to check this. Suffice to say that the asymptotic 
behaviour of (N^) is indeed not enough to find the leading order dependence of the 
BDG action in the flat spacetime interval. 

5 Summary 

We have found a family of causal set boundary terms that agree in the mean with the 
Gibbons-Hawking-York boundary term for a spacelike hypersurface. We presented a 
heuristic argument for how the fluctuations of these boundary terms go with p, and 
provided numerical evidence in 2, 3 and 4 dimensions to support that reasoning. In 4 
dimensions the fluctuations decrease with p. It would be interesting to also find causal 
set analogues of the GHY boundary term for timelike boundaries. The situation is more 
complicated in that case because the identification of “nearest neighbours” to a timelike 
hypersurface in terms of causal structure is less straight-forward than in the spacelike 
case. 

The other major result of this paper is that the average over sprinklings of the 
BDG action for an interval in Minkowski spacetime is proportional to the volume of the 
“joint” of that interval. There is still more work to be done to determine what sort of 
boundary contributions are already contained in the BDG action. 

Other interesting results obtained along the way were causal set expressions for the 
spatial volume of a spacelike hypersurface and the dimension of the manifold the causet 
has been sprinkled into. 
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